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Abstract. We report the first likely spectroscopic confirmation of a z ~ 10.0 galaxy from our ongoing search for 
distant galaxies with ISAAC/VLT. Galaxy candidates at z > 7 are selected from ultra-deep JHKs images in the 
core of gravitational lensing clusters for which deep optical imaging is also available, including HST data. The 
object reported here, found behind Abell 1835, exhibits a faint emission line detected in the J band, leading to 
z — 10.0 when identified as Lya, in excellent agreement with the photometric redshift determination. Redshifts 
z < 7 are very unlikely for various reasons we discuss. The object is located on the critical lines corresponding to 
z=9 to 11. The magnification factor \i ranges from 25 to 100. For this object we estimate SFR(Lya) ~ (0.8 — 2.2) 
M©yr _1 and SFR(UV) ~ (47 — 75) M© yr _1 , both uncorrected for lensing. The steep UV slope indicates a 
young object with negligible dust extinction. SED fits with young low-metallicity stellar population models yield 
(adopting \i — 25) a lensing corrected stellar mass of M* ~ 8 x 10 6 M©, and luminosities of 2 x 10 10 L©, 
corresponding to a dark matter halo of a mass of typically Mtot ^ 5 x 10 8 M©. In general our observations show 
that under excellent conditions and using strong gravitational lensing direct observations of galaxies close to the 
"dark ages" are feasible with ground-based 8- 10m class telescopes. 
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1. Introduction 

Spectacular progress during the last decade has permitted 
direct observations of galaxies and quasars out to redshifts 
z ~ 6.6 (Hu et al. 2002, Fan et al. 2003, Kodaira et al. 
2003), or in other words over more than 90 % of cosmic 
time. With the exception of the cosmic microwave back- 
ground, direct exploration of more distant objects has so 
far been hampered by the need for observations beyond 
the optical domain, and technical difficulties related to 
near-IR observations of such faint objects. However, by 
combining a strong magnification of background galaxies 
by well known gravitational lensing clusters with present- 
day near-IR instruments on 8- 10m class telescopes such 
an attempt appears feasible (Schaerer & Pello 2002; Pello 
et al. 2003). We present here the first spectroscopic results 
obtained from our prototype observational program with 
ISAAC/VLT on one of our best z ^ 7 galaxy candidates. 

In Sect. we summarize the observational strategy 
adopted in the present study, and we present the pho- 
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tometric and spectroscopic characteristics of this object. 
The redshift identification is discussed in Sect. OH The 
physical properties of this galaxy and implications of our 
finding are briefly discussed in Sect. 21 



2. Strategy, observations, and redshift 
determination 

To search for galaxy candidates at very high redshift 
(z > 7) we used the following strategy. We apply the tra- 
ditional Lyman drop-out technique (Steidel et al. 1996) to 
deep optical images. At redshift z > 6 the main spectral 
discontinuity is due to the Gunn-Peterson (GP) trough, 
i.e. the nearly complete absorption of the flux shortward 
of Lya (1216 A), due to the large neutral H column den- 
sity in the intergalactic medium (IGM). In practice we 
search for objects which are non-detected in optical (VRI) 
bands. In addition we require a detection in J, H, and Ks, 
a fairly red J — H due the appearance of the GP trough in 
J, and a blue H — Ks colour corresponding to an intrin- 
sically blue UV restframe spectrum such as expected for 
UV bright starbursts. These criteria, devised from simu- 
lations based on synthesis models at different metallicities 
(Schaerer & Pello 2002, Schaerer 2003, Pello et al. 2003) 
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Fig. 1. Composite JHKs ISAAC image showing: a) The 
core of the lensing cluster A1835, with the position of 
the long slit used during our spectroscopic survey with 
ISAAC/VLT, together with the location of objects #1916 
and #2582, the reference acquisition star (circles), and 
a nearby field galaxy seen on the 2D spectra. Large and 
small dots show the position of the external and internal 
critical lines at z = 10. b) Thumbnail image in the HST- 
WFPC2 R band, where the object is not detected, c, d, e) 
J, H, and Ks images respectively, f) Composite JHKs 
image. All thumbnail images are scaled to the original 
ISAAC pixel size, without smoothing, and are displayed 
with a linear scale. 

allow us to select galaxies with redshifts between 7 and 10. 
Furthermore the detection in at least two bands longward 
of Lya and the combination with the above H — Ks colour 
criterion allows us to avoid contamination by cool stars. 

2.1. Photometry 

Ultradeep JHKs imaging of the central 2x2 arcmin 2 of 
the lensing cluster Abell 1835 (z = 0.253) was obtained 
at ESO/VLT with the near-IR imager and spectrograph 
ISAAC in February 2003. The total exposure times and 
the - excellent - seeing conditions are given in Table 1. 
Indeed, a large fraction of the total exposure time in H and 
Ks was obtained under excellent seeing conditions, with 74 
% of the data below 0.6" in H and 84 % below 0.4" in Ks. 
Photometric data were complemented by deep VRI ob- 




Fig. 2. True color composite images displaying the color 
difference between the lensed galaxy #1916 and other ob- 
jects in this field: a) JHK and b) RJH. This object is blue 
in JHK as compared to cluster galaxies due to its intrin- 
sically blue UV rest frame spectrum in HK, whereas it is 
red in RJH due to the presence of a break within the J 
band. 



servations taken at the Canada-France Hawaii Telescope 
(see Table 1) and R band (F702W) WFPC2/HST images 
(Programme ID:8249, PI: Kneib, exposure time 7.5 ksec). 
Near-IR images were combined and stacked following a 
procedure similar to the one used by Labbe et al. (2003), 
which is well suited for point sources like the ones ex- 
pected here. Some modifications were adopted to improve 
the background subtraction within a cluster of galaxies. 
Photometry was performed on a 1.5 arcsec aperture using 
the SExtractor package (Bertin & Arnouts 1996). Instead 
of adopting the SExtractor error bars, which are based 
on sky noise only, we have computed more realistic back- 
ground noise statistics by randomly sorting 200 blank re- 
gions of 1.5 arcsec aperture on the neighbouring sky area 
~ 10 arcsec around the object. The resulting error esti- 
mates (adopted in Table 1) are larger by 0.1 mag (J,H) 
and 0.2 mag (Ks) than the SExtractor values. The above 
procedure has also been used to derive the random noise 
within 0.6" apertures used to constrain the limiting mag- 
nitudes in the filters where the object is not detected. We 
have also checked that all the optical and near-IR colors 
for cluster galaxies and, in particular the H — Ks versus 
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Table 1. VRIJHKs photometry of #1916. Near-IR mag- 
nitudes are measured within 1.5 arcsec apertures on seeing 
matched images. Limiting magnitudes in the visible bands 
correspond to a detection limit of 1 a within the near-IR 
reference seeing aperture, which is equivalent to a 3 a 
detection on 4 ISAAC pixels. Table entries denote the fil- 
ter name (col. 1), effective wavelength (2), total exposure 
time (3), the image quality /seeing (4), and the observed 
magnitudes in the Vega (5) and AB (6) systems. 



Filter 


A e ff 


^exp 


seeing 


mag 


mag 




H 


[ksec] 


n 


[Vega] 


[AB] 


V 


0.54 


3.75 


0.76 


> 28.5 


> 28.6 


R 


0.66 


5.4 


0.69 


> 28.5 


> 28.7 


I 


0.81 


4.5 


0.78 


> 27.6 


> 28.1 


J 


1.26 


6.48 


0.65 


25.83 ± 1.00 


26.77 ± 1.00 


H 


1.65 


13.86 


0.5 


23.59 ± 0.25 


25.00 ± 0.25 


Ks 


2.16 


18.99 


0.38 


23.64 ± 0.36 


25.51 ± 0.36 



J — H color-color diagram, are fully compatible with the 
expected values for z = 0.25 galaxies. A more detailed re- 
port of these observations and analysis of the photometric 
data will be published elsewhere (Richard et al. in prepa- 
ration). 

Applying the above search criteria has provided 
six z > 7 galaxy candidates in the observed field. 
Among them, the best candidate (called #1916 here- 
after) was retained for subsequent follow-up spectroscopy. 
The coordinates of this object are a(2000) =14:01:00.06, 
(5(2000) =+02:52:44.1. The location and broad-brand 
spectral energy distribution (SED) of #1916 are shown 
in Figs, n an d El photometric data is given in Table 1. 
The radius of the circles in the thumbnail images in Fig.^ 
is similar to the aperture used to compute magnitudes. 
Limiting magnitudes adopted in Table 1 correspond to a 
detection limit of 1 a within the near-IR reference seeing 
aperture (0.6"), which is roughly equivalent to a detection 
limit of 3 a within a region equivalent to 4 ISAAC pixels. 
The source is not detected on the WPC2/HST F702W 
image with R(F702W) > 27.3 mag at 2 a on 4 HST pix- 
els. Figure [3 shows the JHK and RJH true color compos- 
ite images displaying the color difference between galaxy 
#1916 and other objects in this field, in particular cluster 
galaxies. 

The broad band SED of #1916 was used to constrain 
the source redshift using our public photometric redshift 
code Hyperz (Bolzonella et al. 2000) including numer- 
ous empirical and theoretical spectral templates for both 
galaxies and AGNs, and keeping internal extinction as a 
free parameter. The resulting redshift probability distri- 
bution shows a clear maximum at redshift z p hot ~ 9-11. 
This photometric redshift estimate is mainly corroborated 
by the finding of a strong break of > 3.1 to 3.7 AB mags 
between VRI and H, which has a high significance in- 
dependently of the definition used for the limiting magni- 
tudes. The low J-band flux is attributed to the presence of 
the GP trough located within this filter at these redshifts. 
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Fig. 3. Optical to near-IR spectral energy distribution 
of the galaxy #1916. Shown are the broadband photo- 
metric measurements and associated 1 a uncertainties for 
JHKs. Three a upper limits on 4 pixels are given for 
the non detections in VRI. The solid curve shows a good 
(98.7 % probability) model fit with a synthetic spectrum 
of a young unreddened starburst with redshift z = 10.0 
of metallicity Z = 1/5Z , constant star formation, and 
an age of 3 Myr. The crosses indicate the corresponding 
broad band fluxes of the model. The upper axis shows the 
corresponding Lya redshift. 



2.2. Spectroscopy 

To attempt a spectroscopic redshift determination we have 
observed #1916 in the J band with ISAAC. The obser- 
vations were taken between 29 June and 3 July 2003, us- 
ing a 1 arcsec slit width in order to achieve an optimum 
positioning using bright sources. The results presented 
here were obtained again under excellent seeing condi- 
tions ranging between 0.4 and 0.5 arcsec (about 3 pixels 
in the composite frames). The source is not resolved in the 
spatial direction. Thus, we obtained seeing-limited spec- 
tra for this source with an equivalent spectral resolution 
R ~ 5400. 

To search for possible faint emission lines, we have sys- 
tematically explored the J band, where Lya should be lo- 
cated for objects within the 8 < z < 10 redshift interval, 
with priorities set according to the photometric redshift 
probability distribution. For #1916, the following effec- 
tive exposure times were obtained in the different adjacent 
ISAAC bands: A centr ai ~ 1.193 /mi (3.6 ksec), 1.255 /mi 
(6.4 ksec), 1.315 /mi (10.8 ksec) and 1.365 /mi (18.9 ksec). 
Each band has ~ 600 A width. Thus the observed spectral 
interval covers the range from 1.162 to 1.399 jam; in other 
words the "surveyed" interval corresponds to redshifts z ~ 
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Fig. 4. 2D spectra showing the detected emission line 
of #1916, as well the nearby field galaxy and the [O ill] 
A5007 line of the galaxy #2582 (z = 1.68) used as a refer- 
ence to stack the data sets. 2D spectra are sky subtracted 
and show the spectral region around the emission-line at 
1.33745 /im, leading to z = 10.00175 ±0.00045 when iden- 
tified as Lya(1215.67 A). The line is seen on the 2 in- 
dependent overlapping bands at 1.315 and 1.365 /im. No 
smoothing was applied to these spectra. Both the object 
#2582 and the nearby field galaxy are identified in Fig.l. 




1.336 1.338 1.34 

observed wavelength (/xm) 

Fig. 5. a) ID spectrum of #1916, extracted from the 
composite 2D spectrum of the 1.315 and 1.365 jam bands. 
Dotted lines correspond to the position of the main sky 
OH emission lines. The continuum is not detected. The 
dashed line displays the level corresponding to 1 a of the 
background noise. The displayed spectrum is smoothed 
to the instrumental resolution, b) Gaussian fit of the red 
half of the emission line, convolved with the instrumental 
profile (dotted line). Profiles corresponding to rest frame 
60, 100 and 200 km s _1 are also displayed for comparison 
(solid thick lines). 



8.5-10.5 for Lya. With a spectral resolution for the sky 
lines of R = 3100, corresponding to the instrumental 1 arc- 
sec slit width, the fraction of spectral band lost because 
of strong sky emission lines, mainly OH emission-lines, is 
of the order of 30%. Spectroscopic data were reduced us- 
ing IRAF procedures and conforming to the ISAAC Data 
Reduction Guide 1.5 1 , using the same procedure described 
in Richard et al. (2003). 

The deep exposures resulted in the detection of one 
weak emission line at the 4-5 a level with an integrated 
flux of (4.1 ± 0.5) x 10 -18 erg cm -2 s _1 at a wavelength 
of 1.33745 /im (see Fig. 0] and Fig. EJ, which appears on 
both the 1.315 and the 1.365 /im bands. If identified as 
Lya the observed wavelength translates to a redshift of 
z = 10.00175 ±0.00045 2 , in excellent agreement with the 
photometric redshift determination. 

In the 1.315 /im domain, we have also detected 3 emis- 
sion lines on a secondary target (called #2582) included 
in the same slit as #1916 (Richard et al. 2003). This 
galaxy, a faint i?-band dropout at z = 1.68, shows [O ill] 
AA4959,5007 A and H/3. The weakest of them, H/3, has a 
flux of (6.6 ± 1) x 10 -18 erg s _1 cm -2 . This provides an 

1 http: //www.hq. eso . org/ instruments/ isaac/ index .html 

2 Note that the given accuracy reflects only the line mea- 
surement assuming a rest frame A = 1216 A. Given the lack of 
knowledge of the intrinsic line profile this represents a lower 
limit on the uncertainty on the true source redshift. 



independent verification that emission features as faint as 
5 x 10 -18 erg s _1 cm -2 could have been detected in our 
spectra all the way from 1.25 /im to the end of the J band 
(1.399 /im), with some gaps due to sky features. At wave- 
lengths below 1.25 /im, the detection efficiency becomes 
smaller. For an unresolved line of 5 x 10 -18 erg s _1 cm -2 
a S/N ratio between 2 and 3.5 (from blue to red) would 
be expected at A < 1.25 jam. If present, an unresolved line 
of 10 _1T erg s _1 cm -2 would have been detected, with a 
S/N of 3.6 to 5.5, all the way from 1.162 to 1.399 jam. 

We used the good S/N of the [O III] A5007 line of the 
galaxy #2582 at z = 1.68, which is seen in both the 1.315 
and 1.365 /am overlapping bands, to compute the offsets 
required to stack both data sets. The reference star in the 
same slit, together with a standard star, were used to care- 
fully model the telluric absorptions in these spectra. The 
correction applied on the relevant wavelength domain is 
smooth. Fig. ^displays the 2D spectra obtained on the two 
different overlapping and independent bands. The com- 
bined ID extracted spectrum around the emission-line is 
shown in Fig.0 

3. On the redshift of #1916 

How plausible is the observation of Lya from such dis- 
tant an object? Despite the presence of neutral HI in the 
IGM at very high z, all galaxies with z ~ 6 — 6.6 have 
been found through their Lya emission and several expla- 
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nations exist for a partial transmission of Lya emission 
through the ambient ISM and intervening IGM (Haiman 
2002, Santos 2003) which are also applicable to higher 
redshifts. In #1916 the comparison of UV and Lya lu- 
minosities indicates such a possible partial transmission 
(see Sect. 4). The observed line is hardly resolved and 
barely asymmetric (see Fig. EJ). Best fits are obtained for 
rest frame widths below 60 km s _1 using a simple gaus- 
sian to fit the red half of the Lya emission convolved with 
the instrumental profile. Line widths beyond 200 km s _1 
(similar to those of the sample of z ~ 5.7 galaxies of Hu 
et al. 2003) are excluded with a high confidence level, but 
intermediate widths of the order of 100 km s _1 cannot 
be excluded because of the low S/N ratio in the wings 
of the line. In any case, starburst galaxies show in gen- 
eral a wide range of Lya profiles, especially objects which 
are not selected by their line emission, like #1916. In our 
case, a narrower Lya profile could also be related to the 
rather small inferred dark matter halo mass (cf. Sect. 4), 
corresponding to circular velocities of the order of v c ~ 10 
km s _1 . Small or negligible line asymmetries as observed 
here can be explained within simple models by a vari- 
ety of combinations of parameters including the intrinsic 
emission line width, the size of the surrounding H II re- 
gion, and galactic winds (cf. Haiman 2002, Santos 2003). 
We therefore conclude that the detection of Lya emission 
from such a high redshift object and its observed profile 
are not unexpected. 

How secure is the identification of the single line as 
Lya, i.e. the redshift z = 10.0 ? Although other line iden- 
tifications corresponding to lower redshifts cannot be ruled 
out in principle, they appear very unlikely for various rea- 
sons. 

— First, from the constraints of the broad-band SED all 
solutions with z < 7 have a smaller probability (see 
Schaerer et al. 2004 for more details). In fact possible 
low- z (z ~ 2.5) solutions correspond either to young 
bursts with implausibly large extinction (Ay ^ 2-3 
mag) necessary to supress the flux blueward of the 
Balmer break, or to old stellar populations, for which, 
however, no emission lines are expected. For low-z so- 
lutions invoking a Balmer break the single strongest 
emission line often seen in the corresponding wave- 
length domain is [O n] AA3726. 0,3728. 8. However, this 
identification is excluded as the doublet would easily 
be resolved in our observations. Indeed, if the line was 
[O n] A3727 observed at z=2.59, it would appear as a 
double line with the 2 peaks separated by ~ 16 pixels, 
the two profiles overlapping even for moderate velocity 
dispersions (~ 50km/ s). However, this is not observed. 
Other obvious solutions such as Ha or [N n] at z=1.03 
are difficult to reconcile with the photometric SED and 
the absence of other lines in this interval. 
A rather "exotic" solution at z ~ 2.5, albeit of lower 
probability, can be found with an empirical template 
of a metal-poor blue compact galaxy with strong emis- 
sion lines such as SBS 0335-052. In this case the ob- 



served strong H band flux is dominated by [O III] 
AA4959,5007, H/?, and H7 emission lines, whereas the 
spectral range covered by Ks is devoid a strong lines. 
However, to decrease sufficiently the rest-frame UV 
flux between J and VRI an exceptionally large ex- 
tinction of Ay ~ 3.6 is required. In such a case the 
observed emission line should correspond to lines be- 
tween [O 11] A3727 and ~ 4000 A, including e.g. rel- 
atively strong lines of [Ne ill] A3868 or [Ne m]+H7 
A3968 (cf. Izotov & Thuan 1998). Whereas [O 11] A3727 
can be excluded (cf. above), other identifications could 
be possible. However, more than one emission line 
could then possibly have been detected in the observed 
wavelength range (cf. above). Although quite unusual, 
such an explanation has the advantage of being easily 
testable, as several very strong emission lines would be 
expected in the H band. 

For the reasons just discussed redshifts z ^ 7 seem very 
unplausible. For z ~ 7-10 other possible identifications 
of lines sometimes observed in ehigh excitation objects 
and Lyman break galaxies include lines at A ^ 1900 A 
such as C in] A1909 A, He 11 A1640 A, C iv A1550, and 
N v A1240 A. However, also in this case, more than 
one emission line could possibly have been detected in 
the observed wavelength range (cf. above). 

— The second argument in favour of our line identifica- 
tion is the consistency between the well defined pho- 
tometric redshift and the spectroscopic measurement. 

— Finally, the location of our object on top of the criti- 
cal lines at z ~ 9 — 11 (cf. Sect. 4) represents an ad- 
ditional element supporting our redshift identification. 
It is worth noting that this fact was not considered as 
a selection criterion. 

For these reasons we consider for now the identifica- 
tion of the observed line as Lya as the most likely one. 
Obviously future detections of other emission lines and 
possible imaging and morphological constraints on our 
lensing prediction would represent an important verifica- 
tion of our redshift determination. Also additional very 
deep imaging (e.g. at 1 fim) could help improving further 
the photometric redshift of #1916 and of other candidates 
of our sample. 

4. Properties of the z = 10.0 galaxy and 
implications 

Assuming "concordance" cosmological parameters (Q m = 
0.3, Q A = 0.7, and H = 70 km s" 1 Mpc" 1 ) the age of 
the Universe at redshift 10.0 is just ~ 460 Myr after the 
Big Bang. In other words the galaxy we have detected lies 
at a distance corresponding to only ~ 4 % of the current 
age of the Universe. 

According to the lensing model for Abell 1835 (Smith 
et al. 2002), within modelling errors the object #1916 is 
located on top of the critical lines at z ~ 9 — 11. Therefore 
the magnification factor fi is very large, ranging with a 
large uncertainty from 25 to 100, i.e. 3.5 to 5 mag! As in 



6 



R. Pello et al.: ISAAC/VLT observations of a lensed galaxy at z — 10.0 



principle the likelihood to obtain the largest amplification 
is small we will adopt the minimum value of /i = 25. Once 
corrected for lensing, the intrinsic AB magnitude of #1916 
is ^ 28.5 and 29 in the H and Ks bands respectively 
(cf. Table 1). Clearly, without strong gravitational lensing 
and excellent seeing conditions near-IR photometry and 
spectroscopy of such a faint source would be impossible 
with current 8- 10m class telescopes. 

Because of its location close to the critical line the 
object must be multiply imaged by the cluster. However, 
the identification of the counter images is impossible with 
the present data. Deeper observations with the Hubble 
Space Telescope in the near-IR should allow to confirm 
our current interpretation. Redshift ~ 7-10 sources are 
expected to have virial radius of the order of < 1 kpc, 
corresponding to ^O.l", and thus they are unresolved on 
our present images. HST images with diffraction-limited 
resolution should thus allow to confirm the high-z nature 
of this source (e.g. by the detection of strong tangential 
magnification or multiple images), to constrain the lensing 
configuration through morphological considerations, and 
to determine the physical scales involved in star formation 
activities at such redshifts. 

Estimates of the star formation rate (SFR) can be ob- 
tained from the Lya luminosity, the UV restframe flux, 
or SED fits. For a "standard" Salpeter IMF from 1-100 
M and constant star formation at equilibrium we find 
a lower limit of SFR(Lya) - (0.8 - 2.2) Moyr" 1 and 
SFR(UV) ~ (47 - 75) M yr -1 , both uncorrected for 
lensing, adopting published conversion factors (Schaerer 
2003, Kennicutt 1998). The lower SFR derived from Lya 
is thought to reflect in large part the effects of photon 
"destruction" on the blue side of the Lya emission line 
by scattering in the IGM and losses due to other possi- 
ble effects (dust, ISM geometry). SED fits with constant 
star formation rate (e.g. shown in Fig. 03 for Z = 1/5Z 
and 3 Myr age) yield typically SFR - 20 - 120 M yr" 1 
(0.8 to 4.8 M yr _1 after lensing correction) depending on 
the age after the onset of star formation. Lower SFR val- 
ues would be obtained for "non standard" IMFs favouring 
massive stars. 

In general the shape of the overall UV SED does not 
allow one to distinguish age, metallicity, and extinction as 
these parameters alter the UV spectrum in a degenerate 
way. In the present case the situation is somewhat differ- 
ent, as the UV restframe spectrum is very blue when com- 
pared to other known starbursts (Heckman et al. 1998) 
and to model predictions (Leitherer et al. 1999, Schaerer 
2003) - the UV slope /3 measured from the H and Ks- 
band flux is (3 ~ — 3.8 ± 2.1. Comparisons with observed 
UV spectra of galaxies and model predictions also show 
that the observed UV is unlikely to be significantly influ- 
enced by line emission (if present). From this we can quite 
safely conclude that the extinction must be small in this 
z = 10.0 galaxy. 

SED modeling also allows us to estimate the total 
stellar mass and luminosity involved in this star form- 
ing event. Scaling young burst models with metallicities 



Z > 1/50Z© to the observed SED we obtain M,-2x 10 8 
M (or lensing corrected values of M* ~ 8 x 10 6 M ) and 
luminosities of L ~ 4 x 10 11 L (2 x 10 10 L lensing 
corrected) for the above Salpeter IMF. Larger masses are 
obtained for older bursts. The total current stellar mass 
estimated for constant star formation over ~ 3 Myr is in 
good agreement with the above estimate. Note that, as 
star formation is likely to continue for longer time, the to- 
tal mass of stars to be formed during this starburst event 
is likely larger than the "current" stellar mass. 

The above estimates show that #1916 is comparable to 
or somewhat heavier than the most massive old Globular 
Clusters or the most massive super star cluster observed 
in nearby starbursts (Mandushev et al. 1991, Ho et al. 

1996, Mengel et al. 2002). Assuming typical values for the 
baryonic/dark matter content and for the star formation 
efficiency (fifc/fi m x /sf ~ 0.015) this translates to a to- 
tal mass of the dark matter halo of M tot ^ 5 x 10 8 M . 
Such "massive" halos correspond to the collapse of > 2 
a fluctuations at redshift z ~ 10, where, if metals are 
absent or very deficient, star formation is expected to oc- 
cur thanks to cooling by atomic hydrogen (Tegmark et al. 

1997, Barkana & Loeb 2001). "Massive" starbursts like 
the one found here should contribute to the cosmic reion- 
ization; they are possibly even the dominant contributors 
(Ciardi et al. 2003). 

In short, the observed and derived properties of #1916 
agree well with expectations of a young protogalaxy expe- 
riencing a burst of star formation at z = 10. If we are 
witnessing the first star formation event, therefore po- 
tentially the formation of a massive primordial (so-called 
Population III) star cluster, cannot be established or re- 
jected from the present data. Additional spectroscopic ob- 
servations searching for other possible emission line signa- 
tures (e.g. He II A1640, C iv A1548+51 A or other metal 
lines) are necessary to answer this question. In fact a mea- 
surement of such lines, although presumably intrinsically 
fainter than Lya, could well be feasible since the observed 
( "transmitted" ) Lya flux represents only a relatively small 
fraction of the flux emitted from the source (cf. above). 
The expected flux in other lines could therefore become 
comparable to the observed emission line flux. 

Finally, our observations show that under excellent 
conditions and using strong gravitational lensing direct 
observations of galaxies close to the "dark ages" are pos- 
sible with ground-based 8- 10m class telescopes. We are 
looking forward to the exploration of this yet unknown ter- 
ritory from the ground and with the forth-coming James 
Webb Space Telescope. 
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